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Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It's well 
known that many porphyrins have transition metal atoms, and we have explored transition metal 
atoms bonded to those porphyrin-like defects in N-doped carbon nanotubes. The electronic structure 
and transport are analized by means of a combination of density functional theory and recursive 
Green's functions methods. The results determined the Heme B-like defect (an iron atom bonded to 
four nitrogens) as the most stable and with a higher polarization current for a single defect. With 
randomly positioned Heme B-defects in a few hundred nanometers long nanotubes the polarization 
reaches near 100% meaning an effective spin filter. A disorder induced magnetoresistance effect is 
also observed in those long nanotubes, values as high as 20000% are calculated with non-magnectic 
eletrodes. 



Since thei r discovery by Iijima in 1991, carbon 
nanotubes^l (CNT) have become the subject of intense 
research due to their potential for applications^ such as 
in novel electronic devices.™*' Furthermore, in a seminal 
paper by Tsukagoshi et al CNTs entered the realm of 
spintronics, whereby one envision the possibility of using 
the electron spin, instead of its charge, as information 
carrier.^ In that work the authors demonstrated that the 
spin coherence length of polarized electrons injected onto 
CNTs is larger than 300 nm. Thus, carbon nanotube de- 
vices could be used to manipulate spins in a coherent 
manner. 

The prototypical spintronics device uses spin-polarized 
electrons, which are injected from a source into an unpo- 
larized region and analyzed by a polarized drain. Within 
this arrange men t the so-called giant magnetoresistance 
effect (GMR}ZIS manifests itself by altering - via an ex- 
ternal magnetic field - the relative orientations of the 
magnetic moments of the electrodes. From a practical 
point of view this setup usually involves sandwiching dif- 
ferent materials. An alternative to this has been given by 
Kirwan et al. whereby initially unpolarized electrons are 
scattered by magnetic impurities adsorbed on the sur- 
face of a segment of a carbon nanotube. This way, both 
the electrodes as well as the device itself are made of the 
same material. 

An alternative to this has been given by Kirwan et al?" 
whereby initially unpolarized electrons are scattered by 
magnetic impurities adsorbed on the surface of a segment 
of a carbon nanotubeP This way, both the electrodes as 
well as the device itself are made of the same material, 
thus avoiding issues related to surface matching at the 
interface consequently hindering spurious scattering. 

However, one of the issues concerning the use of CNTs 
as spintronics devices is the need to incorporate dopants 
or defects in order to change their electronic transport 
properties. Closed shell species do not interact very 
strongly with the pristine wall of a carbon nanotube™ 
Furthermore, transition metal atoms are more likely to 




FIG. 1: Images of the a) Porphyrin-like defect in a carbon 
nanotube. b) Heme-B molecule, with an iron atom bonded to 
the four nitrogen atoms, c) Segment of a nanotube containing 
an Heme B-like defect. Colour code: grey - carbon; blue - 
nitrogen; red - iron; orange - oxigen 



form clusters when interacting with the pristine wall of 
the nanotube.^ Even linear chains of Fe atoms are more 
energetically favorable than free standing iron atomsP^ 

One possible path to circumvent this problem is to in- 
corporate doping agents during the growth process. In 
that context, carbon nanotubes sythetized in a nitrogen- 
rich atmosphere - the so-called CN X nanotubes - are po- 
tential candidates E2KH1 it has already been demonstrated 
that these nanotubes could be, for example, used as gas 
sensors for a variety of chemical species™ 

The most stable defect in these structures is a 
pyridi ne-like defect consisting of a 4 nitrogen divacancy 
(4ND)!2°mJ T hi s defect (shown in fig. ^) is formed by 
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two vacancies surrounded by four substitutional nitro- 
gen atoms. We have previously exploited the reactivity 
of this defect to attach ammonia molecules and study 
the behavior of the system as a sensor Interestingly, 
this defect is similar to molecules in the porphyrin class, 
in particular, to a molecule known as Heme-B (shown 
in fig. [l]}), which is found, for instance, in hemoglobin 
and myoglobin. This Heme B molecule has an iron atom 
bonded to the site with four nitrogens. Thus it is intu- 
itive to assume that an iron atom - and other transition 
metal (TM) atoms - gets bonded to the 4ND defect of 
the carbon nanotubes. 

The Heme B-like defect has been recently synthesized 
by Lee et, alP^ Their stability was studied by means 
of repeated cyclic voltammograms, and they have not 
observed significant differences after 10 5 cycles, which is 
attributed to the stability of the covalent incorporation 
of the atoms. 

The use of carbon nanotubes - or any other long one 
dimensional system - with defects, however, poses an ad- 
ditional problem; the position of the defects cannot be 
controlled during growth. The result is a device with 
a large number of randomly positioned defects. Thus, 
in order to obtain quantitatively meaningful theoretical 
predictions, one must take into consideration the effects 
of disorder. 

It is generally assumed that disorder has a detri- 
mental effect. Recently, however, we have shown that 
one-dimensional Boron-doped graphene nanoribbons can 
present near-perfect conductance polarization due to 
disorder.^ Such an effect should not depend on the sys- 
tem under consideration provided there was a difference 
in transmission probabilities for majority and minority 
spin cases for a single scatterer. Furthermore, the intro- 
duction of a large number of defects with non-zero mag- 
netic moment leads not only to structural disorder, but 
also to a magnetic one. In fact, the magnetic moments 
of the impurities might be pointing in random directions. 
We can then consider that, in the absence of a magnetic 
field, approximately half of the moments are pointing up 
and half of them are pointing down. A magnetic field 
would tend to align the magnetic moments leading to an 
analogous to the magnetoresistance effect without the 
need to rely on a multilayered material. 

In this work we show that porphyrin-like CNTs can 
exhibit a nearly-perfect polarization and extremely large 
disorder induced magnetoresistance (as high as 20000%) 
driven by disorder. The disorder was simulated using 
CNTs containing a large number of magnetic impuri- 
ties randomly positioned along the tube. For those 
calculations we u sed a combination of density func- 
tional t heory£ il251 coupled to recursive Green's functions 
methodfp™^ 

We have initially performed ab initio calculations 
within Density Functional Theory (DFT^PMH f or a seg- 
ment of a (5,5) CN X nanotube containing a 4ND de- 
fect (figure [la) . As we have said before, the porphyrin 
molecules can have different TM in the 4ND site, as 



shown in figure [T]d. We, thus, performed calculations us- 
ing iron, cobalt, manganese and nickel atoms in the mid- 
dle of the 4ND defect. The final arrangement for the case 
of iron is presented in figure[T];. As one can see, nine irre- 
ducible cells of the pristine system were used to describe 
the region containing the defect. The computational code 
used was SIESTA^ 30 ! 31 ! which uses a linear combination of 
atomic orbitals (LCAO) as basis set. In the particular 
case of this work we have employed a double zeta ba- 
sis set with polarization orbitals. We used the General- 
ized Gradient Approximation (GGA) as parametrized by 
Perdew-Burke-ErnzerhoP21 (PBE) for the exchange cor- 
relation functional. Finally, the atomic coordinates have 
been optimized using a conjugated gradient scheme until 
the forces on atoms were lower than 0.03 eV/A. 

Furthermore, in order to assess whether the transition 
metal atoms in adjacent cells are magnetically coupled^ 
we simulated a supercell with 18 irreducible cells (twice 
the initial size) with two Heme-B-like defects; one case 
with ferromagnetic ordering (ft)i an d another with a 
antiferromagnetic one (ti). The total energy difference 
(E^-E^) is negligible for all transition metal atoms con- 
sidered in this work, so we infer that there's no magnetic 
coupling in our system. 

For the electronic transport calculations, the system - 
following the procedure proposed by Carolli et aP^ - is 
initially divided in three regions namely, the right and 
the left electrodes and a central scattering region. The 
electrodes for our system are taken as semi-infinite rep- 
etitions of the pristine carbon nanotube. In the absence 
of spin-orbit interactions one assumes the two spin fluid 
approximation, whereby one can calculate the electronic 
transport properties of the majority and minority spins 
indepen dentl y of each other. We then use the Landauer- 
BiittikeJ 34 * 35 ! formula to calculate the transmission coef- 
ficients of the system. 



FIG. 2: A long nanotube built from small pieces randomly 
distributed and with different rotations, there are also pristine 
pieces randomly placed. 

In order to access the electronic transport properties of 
a more realistic system with hundreds of nanometers we 
use a co mbinatio n of DFT and recursive Green's function 
methods) 20 1 21 * 261 ^ To do so, we split up a long nanotube 
into small pieces. Each piece is simulated using a sep- 
arate DFT calculation as already described previously, 
and the Hamiltonian and overlap opperators respectively 
H% and S% are stored. For our (5,5) CN X NT we also 
have to consider five different rotations for the position of 
the defect. With those smaller blocks we build up a long 
nanotube, ranging from 20 to 600 nanometers, by ran- 
domly placing the segments with defects together with 
pristine pieces, as shown in fig. [2] One then recursively 
reduce the system to two renormalized electrodes coupled 
via an effective scattering potential that contains all the 
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information about the central region. 

In the low bias limit the differential spin-dependent 
conductance can be calculated by the Landauer-Buttiker 
formula for the currentP^ 
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where f(E',p) is the Fermi distribution function for a 
given temperature. The total conductance is then given 
by the sum of the majority and minority conductances, 
9 = 5 T + 9 l - 

We are interested in two quantities. Firstly, in order 
to quantify the spin filtering eff ect of this device, we cal- 
culate the degree of polarizatio n 1 37 * 38 ! 
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Secondly, as discussed earlier one also needs to take 
into consideration the relative orientations of the mag- 
netic moments. One can analyze the changes in con- 
ductance due to an external magnetic field that tends 
to align the local magnetization of each impurity. The 
magnetic field in our calculations is taken into consid- 
eration only in the alignment of the magnetic moments, 
it has no other effect on the electronic structure of our 
system. Consequently, the value of this disorder induced 
magnetoresistance (MR) is given by 

MR (%) = gl00% " g50% x 100 (3) 
5ioo% 

where gioo% corresponds to the total conductance for 
the case where all of the magnetic moments are pointing 
along the same direction, and <?5o% is the total conduc- 
tance for the case where there is an equal distribution of 
positive and negative magnetic moments. 

Finally, in these long and disordered CN X nanotubes, 
different defect distributions along the CN X give different 
values of conductance. In order to get statistically mean- 
ingful values of conductance we have calculated about 
200 random arrangements for each concentration and 
length of the nanotubes. 

Upon placing the different TM atoms in the 4ND, we 
observe that they strongly bind with a similar final struc- 
ture in all situations. Table |T] shows the binding energies 
and the final magnetization of the system for each one 
of the TM atoms. As can be seen, the binding energies 
are relatively high (the reference is the isolated atom in- 
finitelly separated from the nanotube). One also sees a 
local magnetic moment in all cases, except for the nickel 
atom. A similar behavior was also observed by Shang, et 
dl. In fig. [3] we present the local magnetic moment, 
771 (t^*) — p^ (r) — p^ (f) of the Heme B-like defect. From 
this figure we can notice a highly localized magnetic mo- 
ment in the iron atom. 

In fig. [4] the spin resolved transmission coefficients as 
a function of energy are plotted for a single 4ND defect 
with cobalt, iron, manganese and nickel atoms. We can 




FIG. 3: Local Magnetization m (r) = p T (r) - p l (r) for CN* 
nanotube containing an iron atom (Heme-B-like defect). The 
net magnetic moment of the system is localized in the iron 
atom. 



TABLE I: Binding energy, magnetization and polarization for 
different transition metal atoms bonded to the 4ND defect. 



Atom 


Fe 


Co 


Mn 


Ni 


Binding Energy (eV) 


-6.3 


-4.5 


-4.2 


-4.2 


Magnetic Moment (pb) 


2 


3 


5 





Polarization (%) 


8.95 


1.46 


2.70 


0.00 



note on fig. [4] different transmission probabilities for the 
up and down spin channels at the Fermi level except 
for nickel, as expected, since nickel shows no local magne- 
tization for this system. For cobalt, iron and manganese 
the electrons with different spins are not be equally scat- 
tered, and this may lead to a spin filter, a system in which 
a non polarized current enters in one side and after pass- 
ing through the system the current is spin polarized. We 
have calculated (with eq. [2]) the polarizations shown in 
table |T] using the conductances obtained for a low tem- 
perature calculation. In all cases there was only a single 
TM defect and one can clearly see that the polarization 
is always bellow 10%. Although the magnetic moments 
are generally high - except for nickel - the localized d 
states are positioned either above or bellow the Fermi 
level. This can seen from fig. [5] where the total density 
of states are shown. In the particular case of iron (figure 
[5})) the minority spin states are closer to the Fermi level 
which consequently lead to a larger polarization.^ Still, 
this leads to values of the polarization which are not very 
high for any of these systems to behave as spin filters. 

For the long disordered CN X we have chosen the Heme- 
B defects (Fe), since they show the highest polarization 
for a single defect. In all cases considered here we have 
taken a defect concentration of 0.65% per mass. In figure 
|6^a-b) we show the spin-resolved logarithm of the conduc- 
tance averaged over more then 200 random realizations 
as a function of device length. As we are in the Ander- 
son localization regime the conductance is supposed to 
vary exponentially as function of length, following the 
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FIG. 4: Transmission coefficients as function of energy for a 
single 4ND defect with different transition metal atoms, a) 
cobalt , b) iron, c) manganese, and d) nickel. 
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FIG. 5: Total density of states as a function of energy for a 
single 4ND defect with different transition metal atoms, a) 
iron , b) cobalt, c) manganese, and d) nickel. 



well know exponential relation 



9 oc e V , 



(4) 



where £ CT is the spin-dependent localization length. This 
can be clearly seen for both the cases of majority and 
minority spin channels. The localization lengths obtained 
by the linear fit are shown in table JTTJ 

As one can clearly see from figurepp the conductance is 
completely spin polarized for the case where all the mag- 
netic moments are pointing in the same direction and 
the nanotube is longer than 200 nm. For the case where 
half the magnetic moments are pointing up and half are 
pointing down the average magnetization is, as expected, 
zero (and therefore not shown in the graph). As shown 
in a previous workJ22J this disorder-driven polarization ef- 
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FIG. 6: Average logarithm for a) majority and b) minority 
spin conductances as a function of length, c) Polarization 
for the 100 % case, and d) average magnetoresistance. All 
graphs shown correspond to a defect concentration of 0.65% 
and temperature of 3K. 



feet is independent of the type of system being studied 
provided there are different scattering probabilities for 
majority and minority spin channels in the single impu- 
rity case. For the 100% case, since £} is larger than 
the spin down conductance will decay much faster than 
the spin up conductance, following 



5 T oc 



(5) 



where n — ^ For the case shown here n = 1.642, 
whereas for the single impurity case, R^/Fft = 1.634. 
These two results are very close, corroborating our pre- 
vious assumptionPSl 

In fig. [6]i we present the values for the magnetoresis- 
tance (calculated with eq[3]as a function of length. The 
calculated MRs in this case reached values up to 20000 %. 
As in the case of the polarization, the presence of ran- 
domly distributed scatterers lead to an enhancement of 
the MR effect up to extremely large values. 

In order to address the effect of a magnetic field applied 
to the system and the possibility that not all the mag- 
netic moments are completely aligned we also considered 
that only 80% of the defects are magnetically aligned. 
The results are shown in figure [7J and the linear fits for 
the localization lengths also are presented in table [TT1 
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Length (nm) 

FIG. 7: Average logarithm for a) majority and b) minority 
spin conductances as a function of length, c) Polarization 
for the 80 % case, and d) average magnetoresistance. All 
graphs shown correspond to a defect concentration of 0.65% 
and temperature of 3K. 

TABLE II: Values of the spin-dependent localization lengths 
£ CT for 0.65% defect concentration and temperature of 3K. 





Localization length (A) 






100% 80% 


50% 


Spin up 


417 357 


302 


Spin down 


254 277 


306 



We can note a decrease (increase) in localization length 



for majority (minority) spin compared to the 100 % case. 
This is to be expected since one is moving towards higher 
magnetic disorder, i.e. the 50 % case. From fig.[7[: we can 
see a polarization near 100% for 700nm long CN X nan- 
otubes. Most importantly, the magnetoresistance (Fig. 
[7ji) presents values which are one order of magnitude 
lower than the fully aligned arrangement, but it is still in 
the 1000 % range. Thus, even in the case where not all 
spins are aligned, there is still an extremely large disorder 
induced magnetoresistance. Thus, this disorder-driven 
GMR effect is extremely robust toward fluctuations in 
the alignment of the magnetic moments. 

We have used an ideal paramagnet model to estimate 
the needed magnetic field to obtain a 80% magnetiza- 
tion at 300K and 3K. Unfortunately the needed magnetic 
field at ambient temperature is about 200 T, making it 
impracticable for ambient temperature devices. For a 
temperature of 3K the needed field will be about 2T for 
an 80% magnetization and about 5T for a 95% magneti- 
zation, so it's possible to observe the predicted effects in 
this paper in low temperature experiments. 

In summary we have observed that transition metal 
atoms bind strongly to nitrogen defects in CN^, carbon 
nanotubes in a fashion similar to Heme B molecules. The 
end result is a scattering site with a localized magnetic 
moment on the transition metal atom that leads to a 
small conductance polarization in the case of a single 
impurity. For a large number of such impurities ran- 
domly distributed along carbon nanotubes a few hundred 
nanometers long, it leads to near perfect polarization and 
a large magnetoresistance (up to 20000%). An interest- 
ing feature of the system proposed here is the fact that 
they do not need polarized electrodes as it is usually the 
case. We estimate that, at low temperature, this effect 
could be measured experimentally. 
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